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Abstract: Since network motifs are an important property of networks and some 
networks have the behaviors of rewiring or reducing or adding edges between old vertices 
before new vertices entering the networks, we construct our non-randomized model N{t) 
and randomized model N'{t) that have the predicated fixed subgraphs like motifs and 
satisfy both properties of growth and preferential attachment by means of the recursive 
algorithm from the lower levels of the so-called bound growing network models. To 
show the scale-free property of the randomized model we design a new method, 

called edge-cumulative distribution, and democrat two edge-cumulative distributions of 
N{t) and N'{t) are equivalent to each other. 
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1 Introduction and concepts 

Barabasi and Albert [3] observed that both ingredients of growth and preferential attaehment are 
needed for the development of the stationary power-law distribution, since growth and preferential 
attachment are mechanisms common to a number of complex systems, including business net¬ 
works, social networks (describing individuals or organizations), transportation networks, and so 
on. Directed networks have been discussed in [3], [8], |16] and |17] . The authors in [5] show ver¬ 
tex replacement and vertex addition methods to produce small-world and scale-free networks from 
a low diameter “backbone” network; specific recursive scale-free constructions with fixed degree 
distributions were investigated in [I2] and [6]. 

Network motifs are an important local property of networks, and have been identified in a wide 
range of networks across many scientific disciplines and are suggested to be the basic building blocks 
of most complex networks BM- Motifs are of notable importance largely because they may reflect 
functional properties, and they have recently gathered much attention as a useful concept to uncover 
structural design principles of complex networks, and may provide a deep insight into the network’s 
functional abilities m- If the probability of a given subgraph to appear at least the same number 
of times as in the real network is smaller than a given threshold, the subgraph is considered a motif 
of the network [7] . In |10l [T^ [TT] the authors defined “network motifs” patterns of interconnections 
occurring in complex networks at numbers that are significantly higher than those in randomized 
networks, and they found such motifs in networks from biochemistry, neurobiology, ecology, and 
engineering. 

We are motivated from the Barabasi-Albert (BA) model that is an algorithm for generating 
random scale-free networks using a preferential attachment mechanism, and motivated from a 
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phenomenon that some groups of new vertices enter simultaneously into a network, not one by 
one. So, we will show recursive construction methods to generate models that are graphs having no 
multiple edges between the same pair of vertices and no self-edges that connects the same vertex 
at both ends, and some edges of the models can be removed randomly, and then some new edges 
can be added to the remainder randomly. For the purpose of simulation, our deterministic models 
are visual and satisfy both properties of growth and preferential attachment, and furthermore 
the construction of our models can be shown as the recursive scale-free algorithm from the lower 
levels of the so-called bound growing network models. Our models differ from those models in 
P [la [ig [H [201EU, since we grow them by adding some graphs called seeds and by making 
edges removed/added randomly. We use these seeds as motifs appeared in many hierarchical and 
biological networks. But these seeds are not directed in this paper. 

As reported, some networks have the behaviors of rewiring or reducing edges between old vertices 
before new vertices entering the networks. Based on this observation, we have developed our 
evolutionary models that we will investigate in this article: one is called the uniformly (r,F)- 
growing network models (abbreviation as (r, F)-ugnms) that have no edges removed/added; the 
another one is called the randomized e-bound growing network model (abbreviation as randomized 
EBGN-model). Here, our “e-bound growing” is the same growing way as that shown in [21], and is 
local not universal as p. Our first step is to generate randomized EBGN-models. The generative 
method for obtaining a high level randomized EBGN-model is based on two phrases: (1) add 
new graphs like motifs to the so-called bound-edges of a lower level randomized EBGN-model; (2) 
the randomized mechanism: based on a probability pr remove some edges from this lower level 
randomized EBGN-model, and then add new edges to the remainder by another probability pa- 
Our second step is to show the scale-free property of the randomized EBGN-models by means of 
the topological structure of (r, E)-ugnms. The term “removing edges” mentioned in this paper 
means that we delete edges and keep their ends in networks; and the term “adding edges” means 
that we join two vertices in networks by an edge if there is no edge between them. Glearly, it is not 
easy to obtain the degree spectrum of a randomized EBGN-model. As an alternative approach, 
we defined the edge-cumulative distribution for proving the scale-free formation of the randomized 
EBGN-models. 

2 Models having edges removed/added randomly 

Our randomized EBGN-model N'{t) having edges removed/added will be constructed by the fol¬ 
lowing deterministic algorithm-I. For short writing we use “i € [m, n]” rather than “i = 

1,... ,n”, where integers n > m > 1; the number of elements of a set X is denoted as |A|; the 
length of a path P is denoted as |P| in the whole paper. An edge uv has its own ends u and v. 
Let F be a seed set of finite connected graphs having (> 1) vertices and me (> 0) edges, and 
let every connected graph G £ F have its vertices xi,X 2 , - ■ ■ ,Xm^. Here, the graphs in F play the 
role of motifs or communities appeared in many real networks. We take an integer r that holds 
> r > 1, and call r the bound thickness as wall as every G £ F a seed. Let N'^^ and 

^be(t) be the numbers of vertices, edges and bound-edges of the randomized EBGN-model N'{t), 
respectively. 

Algorithm-I 

Step 1. (Initialization) The initial model A'(0) has no multiple edges and loops, and 
is connected. A'(0) has its own vertex-set F'(0) with N^^ = |F'(0)| > 2 vertices 
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and its own edge-set E'{0) with iV'g = |i?'(0)| > 1 edges. We define every edge 
uv of N'{0) as a bound-edge. Let B'{0) be the bound-edge set of N'{Q). Clearly, 
B'{R) = £'^(0). For t = 1, the new model A^'(l) can be obtained by doing the 
following operation-I. 

Operation-I: Add a seed G a F to each bound-edge uv of A^^(O), and join every vertex 
Xi of G with vertex u and vertex v to form two edges XiU and XiV, respectively; 
and select randomly vertices , Xi 2 , ■ ■ ■, Xi^ to define 2r edges Xi. u and vxi. as the 
bound-edges for j € [l,r]. 

Step 2. (Iteration including growth and removing/adding edges) For t > 2, N'{t) is 
obtained by doing the operation-I to each bound-edge of N'{t — 1); remove randomly 
some edges from — 1) by the probability pr, and the remainder is denoted as 
T'(t — 1); and then add randomly new edges to T'{t — 1) according to the probability 
Pa, where 0 < pr,Pa < I and pr is independent of pa- 


Notice that N'{t — 1) is not a subgraph of N'{t) by the construction of the Algorithm-I. We call 
the Algorithm-I the growth-first randomness-second algorithm. An illustration about the model 
N'{t) is shown in FiglU Figj2]and FigjSj We have the following basic parameters 


^'v,t — + '’^vN'a Q ^(2r)^ — 0 + ^ _ 1 ’ 

k=0 


^e,t — P* + 0 

NM = {2rYN'. 


, (2r)* — p* 


2r — p ’ 


( 1 ) 


wherep = (I—pr)(l-|-pa)- Here, two numbers N^^^ and Nfi^ are shown in (l5|), and the edge number 
Ng j will be deduced in the proof of Theorem [1] in Appendix. 



Figure 1: (a) Two seeds with = 4 and mg = 3; (b) iV'(O); (c) iV'(l) having 15 vertices and 36 edges; 
(d) the procedure of growing iV'(l) having 63 vertices and 188 edges by the bound thickness r = 2. 
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(e) (f) 

Figure 2: (e) T'(l) is obtained by removing 18 bold-edges of N/1) shown in FiglUby the probability 

(f) the randomized EBGN-model N'{2) obtained by adding 9 new dashing-edges to T'{1) by the probability 

Pa = ^- Here, the number of edges of A^'(2) is less than that of N'{1). 


Let {k)Q = be the average degree of the initial model N'{0). By Eq.([I]), the average degree 
{k)' of N'{t) is easily estimated as follows. 


{k/ 


(me-b 2m„)(2r - 1) 

- — (X - 

N/^ my{2r-p) 



2r-l 
2r — p' 


t —>■ oo. 


( 2 ) 


The result (l2|) shows that N'{t) is a tree-like model. Since pr = number of removing edges 

number of adding edges is equal to pa{I — Pr)N/^_i, we announce that pa > as 0 < pa < 
and 0 < Pa < when ^ < p^ < 1 for the existence of the randomized EBGN-model N'(t). We 
define the edge-cumulative distribution P/cumi^) the model N'{t) by 


K-c.Jk) = /r(VKA 

with 0 < r < t. We can show a result in Theorem [1] whose proof is detailed in Appendix: 

Theorem 1. The randomized EBGN-model N'(t) is scale-free, and its edge-cumulative distribution 
P'e-cum^^) obeys the power law distribution and is independent of p = (1 — pr)(l +Pa)- 


3 Models having no edges removed/added 

We construct our (r, T)-ugnms N(t) for time steps t > 0 in the following deterministic algorithm-II. 
The seed set F is defined well and the bound thickness r subject to m„ > r > 1 in Section 2. 

Algorithm-II 

Step 1. (Initialization) Eor t = 0, A(0) is a connected network model having no 
multiple edges and loops. A'(O) has its own vertex-set E(0) with N^^ = |F(0)| (> 2) 
vertices and edge-set E{0) with Nefi = |E(0)| (> 1) edges. We define every edge 
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(g) 


(h) 


Figure 3: (g) T'(l) is obtained by removing 12 bold-edges of N/l) shown in FiglT] shown in FiglT]by the 
probability Pr = ^; (h) the randomized EBGN-model N'{2) obtained by adding 16 new dashing-edges to 
T'(l) by the probability pa = §• In this case, the number of edges of N'{2) is greater than that of N'{1). 


of A^(0) as a bound-edge and call A^(0) the initial {r, F)-ugnm. Let B{0) be the 
bound-edge set of N{0). 

Step 1. (Iteration) For t > 1, an (r, F)-ugnm N{t) is obtained by doing the following 
operation-II: 

Operation-II: Add a seed G G F to each bound-edge uv of N{t— 1), and joining every 
vertex Xi of G with vertex u and vertex v, respectively, produces two edges XiU and 
XiV\ and select arbitrarily vertices xi^^xi^, ... ,Xi^ to define 2r edges Xi-u and vxi^ 
for j G [l,r] as the bound-edges of N{t). 

Thereby, we write X{t) = V(t) \ — 1) to indicate the set of new vertices added into N{t — 1), 

Y{t) = E{t) \ E{t — 1) the set of new edges added into N(t — 1), and B{t) the set of bound-edges 
of N{t). Clearly, different bound-edges of N{t — 1) may correspond to different seeds in F. It is 
not difficult to observe that N{t — 1) is a sub-model of N{t) for t > 1, so we also call N{t) a nested 
(r, F)-ugnm. 


3.1 Basic properties 

We denote the degree of a vertex x of an (r, F)-ugnm N{t) at time step t as deg(a;,t) that is the 
number of edges incident to x in N{t), and write the numbers of vertices, edges and bound-edges 
of N{t) by Wy = |F(t)|, Ne,t = \E{t)\ and Nbe,t = respectively. 

Basic numbers. By the construction of an (r, F)-ugnm N{t), we are not difficult to obtain 


W,1 = N^,0 + myNefi, W,1 = Nefl + {me + 2my)Nefl, Nbe,l = 2rNbe,o, 


where Nbe,o = A"e,o- For k > 2, we have the following recursive formulas 


(3) 


^v,k — ^v,k—l “t“ 

^e,k ^e,k—l (^e “1“ 
^be,k — ‘^'^^be,k—l‘ 


( 4 ) 
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Furthermore, the numbers N^^t, N^^t, N^e^t of vertices, edges and bound-edges of an (r, F)-ugnm 
N{t) at time step t > 1 are as follows: 

(2rY — 1 

Ny t — Ny Q TTlyN^ Q - - , 

2r — 1 

Ne,t = Nefl + (me -h 2m^)iVe,0^^-^-7^, 

2r — 1 

Nhe,t = (2r)*iVe,o. 


Eq-Q shows N{t) has the growth behavior as t tends to infinity. 

Average degree. The average degree (k) of an (r, F’)-ugnm N(t) at time step t > 1 can be 
estimated as 


{k) 


“^Ne^t _ 2 ^e,0 + (me + 2m^)Nefl Z]fc=o(2^)^ 

~ ■ A,,o + m,Ae,oEfc=o(20" 


oc 


2(me + 2m„) 
rriy 


4 + {k/o, 


( 6 ) 


for larger t, where (A:)o = average degree of the initial (r, F’)-ugnm A^(0). Clearly, an 

(r, F’)-ugnm N[t) is a sparse model for larger t, since {k)o is a constant. Also, N{t) is tree-like, 
since the number N^^t of edges of N{t) grows linearly with the number of vertices of N{t) at 
time step t. 

Notice that two numbers of the newly added vertices and the newly added edges in forming an 
(r, F’)-ugnm N(t) are 


\X{t)\ = \V{t) \V{t-l)\= m,(2r)*-iAe,o, 

|y(t)| = \E{t) \ E{t - 1)1 = 2|X(t)| = (me + 2m,)(2r)*-iAe,o, 


respectively. The new increasing part of an (r, T)-ugnm N(t) at time step t is tree-like, that is, the 
number |T(t)| of new edges added to N(t — 1) grows linearly with the number |A(t)| of vertices 
newly added according to oc 2 as t tends to infinity. 


3.2 Degree spectrum, diameter 

By the linear preferential attachment rule in Ref. [7] , the vertices having larger degrees in N (t) play 
the role of connecting network and attract easily new vertices that are entering into N{t). A vertex 
w added newly is easily adjacent to those vertices z of the network N{t) having larger degrees, that 
is, 

deg( 2 ;,t) deg{z,t) 


Pt{w z) = 


2N, 


e.t 


( 8 ) 


EveV{t) degiv,t) 

So, finding the degrees of vertices is very important for figuring out many properties of networks. 
For a vertex u in the initial model A(0), it is not hard to compute its degree in N{t) 


1 + m^ 


H - 1 


r — 1 
(1 -|-tm,;)deg(ti,0) 


deg(«,0), r > 2; 

r = 1. 


deg{u,t) = 

[ (1 

And the vertex u is the common end of r*deg(u, 0) bound-edges of A(t) when t > 1. 


(9) 
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In forming N{s) at time step s > 1 with respect to s < t, a vertex a; of a seed G € F is added 
newly into N{s — 1) and is selected as an end of some bound-edge of N{s). So, in N{t) at time 
step t > s > 1, the vertex x has its own degree 


deg(x,t) 


r* ^ — 1 

deg(3:, G)-I-2-I-2m„-—, r > 2; 

r — 1 

deg{x,G)+ 2+ 2mv{t - s), r = 1, 


( 10 ) 


and is the common end of 2r*“^ bound-edges of N(t). 

For a vertex y of a seed G £ F being not an end of any bound-edge of N{t), we know its own 
degree deg(y, t) = deg(y, G) + 2. Eq.s ([8]), Q and (fT0|) show that our (r, F)-ugnm N{t) possess the 
behavior of preferential attachment, since for larger t and r > 2, we have 


Pt{w u) > Pt{w —> x) > Pt{w —> y). 


We can rewrite new added-vertex set X{s) at each time step s € [l,t — 1] with t > 2 by 
X(s) = X’^/ U Xf/ where X’^/ is the set of newly added vertices that are selected as the ends 
of some bound-edges of N{s), and every vertex of the set = W(s) \ X^^ is not an end of 

any bound-edge. It is not difficult to compute both numbers \X\'^\ and by Eq.(I7|). We 

can see that \Xl^\ = rA^e,o(2r)^“^, and the degree of each j/j of X^^ in N{t) is shown in Ea. (fT0]l : 
\X/^^\ = {rriy — r)N(,fl{2rY~^, each yj of X/^^ has its own degree deg{yj, G) -|- 2 in N{t). 

By the above degree spectrum of N{t) we can compute the probability P{k) that a randomly 
selected node has exactly k edges. Suppose that each vertex of X^^ has its own degrees greater 
than k for s G [1,t] with 0 < r < t, and P{k* > k) is the probability of vertices having degree not 
less than k. So, by Q and (fTO)) . we have deg(u, t) > k for u € V (0), and deg(x, s) > k for x G 
with s G [l,r], but X^^ contains no vertex x having its own degree deg(x,s) > k. Thereby, 


P{k* >k) = 


1 


K,t 

Xv,0 + fXefl 


’'be I 

*■5 I 


S=1 

{2rY - 1 


1 


2r - 1 


N,flPrNe,oY.{2ry-^ 

^f=i 

(2r)* - 1 


N.,t 

Xy 0 “ 1 “ N^ nTTly 


2r - 1 


oc —(2^)"-^ 
rriy 

and P{k* > k — 1) (x ^(2r)'^+^“*, which enable us to get 


P{k) = P{k* >k-l)- P{k* >k)(x ^ - ^{2ry-\ 


r{2r — 1) 


rriy 


Let h{r) = Plugging t = t- h{r) Ink into Eq.([l2]) can show 

p^k) oc !:^LzJlk-Hr)H 2 r)^ 


( 11 ) 


( 12 ) 


since 


{2rY~^ = (2r)“^l^P“^ = . 


Hence, N{t) is an exponentially growing model, also, is a scale-free model as mentioned in [T]. 


( 13 ) 
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Diameter. The notation dis(x, y) indicates the distance between two vertices x and y of 
N{t), and D{t) denotes the diameter of N{t) at time step t. Note that dis(x,y) is the length 
of a shortest path connecting x and y, and D{t) = max{dis(x,y) : x,y ^(^)}- Foi' ^ iiew 
added vertex x £ X{s) at time step s > 2, the algorithmic construction of N{t) tells us that x is 
adjacent to two vertices Ui and Uj, where Ui £ X{i) and Uj G X{j) with i j such that {i,j) £ 
{(0, s — 1), (1, s — 1 ), ( 2 , s — 1 ),..., (s — 2 , s — 1 )}, and X(0) = 1^(0). So, x can be connected with a 
vertex uq of the initial model N ( 0 ) by a path xUs- 2 Us-a ■ ■ ■ U 2 Uq or another path xUs-iUss • • • uiuo- 
Thereby, two vertices x' and x" of N{t) can be connected with two vertices uq and vq of A^(0) by 
two paths P(x',uo} = x'ui- 2 Ui -4 ■ ■ ■ U 2 U 0 and Q(x",vo} = x"uj- 2 Uj -4 ■ ■ ■ U 2 V 0 , respectively. Let 
U(uo,vo} be the shortest path connecting uq and vq in iV(0), so its length |C/(tto,no)| < D(0). We 
can estimate the distance dis(x',x'^) in the following inequalities 


dis(x',x") < \P{x',uo)\ + |C/(tto,no)| + \Q{x'',vo)\ < 2 


+ ^( 0 ), 


this result shows D{t) < t + 1 + D(0) at time step t. Therefore, N(t) is small-world because its 
diameter D{t) is the same rank as logNy^t, that is, D{t) = 0(log W,t)- 


3.3 New properties 


We present some new statistical approaches for exploring connections between known and new 
statistics methods. 

Newly added average degree. We can get the newly added average degree {k)new as follows. 


_ 2|y(t)| _ 2|£;(t)\j^(t-i)| 

\x{t)\ |y(t)\y(t-i)| 


2{me + 2my){2rY ^jVe,o 
my{2ry-^Nefl 


4 {k)o oc {k). 


(14) 


Clearly, two averages {k)new and {k) are the part and the whole, respectively; so that the part 
{k)new is equivalent to the whole (k) according to Eq.® and Eq. lfTTll for sufficiently large t. 

Edge-cumulative distribution. Eor an integer 6 with respect to 0 < d < t, we define the 
edge-cumulative distribution T’e-cum(^) of an (r, F)-ugnm N{t) by Pp- mm (k) = 7 ^ X]f=oBy 
Eq.([S]) we have the following result. 


Theorem 2. The edge-cumulative distribution Pe-cum{k) of the (r, F)-ugnm N(t) obeys the power 
law distribution, since 

2 r — 1 

The proof of Theorem [2] is developed in Appendix. 

When A^(0)’s vertex number my > 1 and the bound thickness r = 1 we obtain 7 ^ = 2-|- ) 


with 2 < 7 fc < 3. Clearly, the case r > 2 shows h{r) ln(2r) = 2r -|- in(' 2 r+m ) ^ other words, 

the edge-cumulative distribution Pe-cum(^) will be influenced only by the parameter r. We claim 
that there is a function /(r, m„) such that P{k) = f{r,my)Pe-cnm{k), according to the following 


Pik) 


Pe-c\im{k) 


OC 


( 2 r - 1/ 

2m„ 


or P{k) oc 


( 2 r - 1 )^ 
2 m,, 


/k). 


(16) 


Observe that P{k* < k) and Pe-cum(^), where P{k* < k) is the probability of vertices having 
degrees less than k + 1. By P{k* < k) = 1 — P{k* > k) oc 1 — 7 ^( 2 r)'^“*, we have 

P{k* < k) (my — r)(2r — 1) 
oc — 


p 

P.-C. 


i(^) 


2 rm„ 


t ^ 00 . 


(17) 
















4 DERIVATIVE MODELS 


9 


The (ffc, efc)-models. Based on the degree spectrum of an (r, F)-ugnm N(t}, we will compute 
the number Sjv(< k) of vertices having degrees no more than k in an (r, T)-ugnm N{t) and the 
sum Qn{^ k) of degrees of these vertices at time step t. Clearly, the number of vertices having 
degrees greater than k is equal to Sn{> k) = k) and the number of degrees of the 

vertices whose degrees are greater than k is equal to Qn{> k) = 2N(,^t — Qn{< k). Here, we 
consider a particular selection as forming N{t). Suppose that every graph G £ F has the vertices 
yi,y 2 , ■ ■ ■, yrriy and a vertex yi has its own degree deg{yi,G) = di for i G [1, m^] such that dj > dj+i 
for i £ [j,my — 1]. In the procedure of constructing an (r, F)-ugnm N{t) by adding a seed G to 
each bound-edge uv oi N{t — 1), we select the previous vertices yi,y 2 , ■ ■ ■ ,yr of the seed G and 
dehne 2r edges yju and vyj as the bound-edges of N(t) for j £ [l,r]. 


Theorem 3. For r >2, two quantities Vk = and e^ = ^ 2 N^t'^ the power law distribu¬ 

tion according to 


Vk oc ^k-^{r)ln{2r) ^ 
rriy 


Ck oc 


rm,. 


m* -\- 2r 


uie + 2m^ (r — l)(2r — 1) 


^-h(r)ln(2r) 


nif, + 2my r — 1 


where h{r) = ^ and m* = deg(yi, G). 

The proof of Theorem [3] is presented in Appendix. 


4 Derivative models 

Our uniformly (r, T)-growing network model N{t) can be used to build up another class of ran¬ 
domized model M'{t) having partially rewiring edges in the way that we do an operation-II to each 
bound-edge of M'{t — 1) and next, rewire some edges of M'{t — 1). Notice that the operation-II 
guarantees the growth and preferential attachment of both models N{t) and M'{t), and two mod¬ 
els both have the same number of edges. By means of the edge-cumulative distribution, we can 
demonstrate that M'[t) obeys the power-law distribution. 

One more randomized bound-growing network model M"[t) can be dehned by applying the 
step 1 of the Algorithm-I, and then, for t >2, M"{t) is obtained by doing an operation-I to each 
bound-edge of M"{t — 1) and make the number M"^ of edges of M"{t) holds V'^ (or 

FIe,t < — ^et)- scale-free behavior of M"{t) can be proven through comparing three 

edge-cumulative distributions of M"{t), N'(t) and N{t). 

5 Conclusion 

We, by means of the topological structure of our uniformly (r, T)-growing network model N(t), 
show the scale-free behavior of our randomized EBGN-model N'{t). However, it is not easy to 
obtain the exact values of some distributions of the randomized EBGN-model N'{t), such as the 
power-law degree distribution, the clustering coefficient distribution and diameter, since we are not 
able to calculate its degree spectrum and N'{t) may be disconnected. To prove the randomized 
EBGN-model N'{t) to be scale-free, we use the principle of comparison to show that two edge- 
cumulative distributions Pe-cum{k) and Pe.cumi^) both are equivalent to each other. In contrast, 
the removing-edge probability pr and the adding-edge probability pa hold pr < ^ < Pa < I, the 
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number of edges of the randomized EBGN-model N'{t) is greater than the number of edges of 
the uniformly (r, F)-growing network model N(t), that is, N'^^> and moreover 

when Pa < ^ < Pr < ^- Likewise, we observe that N'{t) has a giant component that charges the 
growth and preferential attachment and makes N'{t) to be scale-free; the other part of called 

the subpart, does not determine the topological structure of N'{t) no matter removing edges from 
the subpart or adding edges to the subpart. We use a notation G{t — 1) to denote the remainder 
after removing and adding edges to N'{t — 1) in the Algorithm-I. It may occur a phenomenon that 
G{t — 1) is regular rather than scale-free at some time step t > 1, that is, G{t — 1) is a random 
Erdos-Renyi type graph that follows the form of the Poisson variable. 

We point out the fragility of our randomized EBGN-model N'{t). If we remove the vertices of the 
initial model A^'(O) from the remainder may have two or more components that are connected 

to each other. For example, attacking the vertices of the initial model A^^(O) at time step to may 
make the remainder M{to) of N'{tQ) has some components for i = 1,2, ...,m (> 2). 

However, according to our construction algorithm of N'{t), a component Mi(to;t) having some 
bound-edges can be self-growing and obeys the power-law distribution for enough larger t > to, so 
that the model M(to; i) based on the initial model M{tQ) will become a scale-free model as adding 
edges make the model M{to] t) to be connected. 

Glearly, the models mentioned above can be generalized into directed models. As further work 
we propose: 

(1) If, to obtain we do firstly removing/adding edges to N'{t — 1), and then grow the re¬ 

mainder by doing the operation-I to each bound-edge of the remainder. This construction algorithm 
is refereed as randomness-first growth-second. Is the model N'{t) scale-free? 

(2) In the step 2 of the Algorithm-I, we substitute by pr^t and pa^t two probabilities pr and pa at 
time step t, such that there are pr,i 7 ^ Pr,j or pa^s 7 ^ Pa,i for some i 7 ^ j or s 7 ^ /. What topological 
structure does N'{t) have? In other words, we want to optimize randomized EBGN-model N'{t), 
as proposed in [ 2 ]. 
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Appendix 


Proof of Theorem [2l To calculate the edge-cumulative distribution of N{t), we take <5 with 
respect to 0 < <5 < t for larger t. So 
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Plugging 6 = t- h{r) Ink for h{r) = + i„( 2 r+m„) Eq.([l8D yields 


P (k) (X h-Hr)H2r) 


2r - 1 

by Eq. (fT^ . as desired. 

Proof of Theorem [3l For r > 2, we can estimate 
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Let m* = deg(yj, G). We can estimate 
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Notice that 2 t(m”+ 2 m ) 0 as t oo. Plugging t = t — h{r) Ink into Eq.([22]) produces 




Qjv(> k) ^ , j^-h(r) ln(2r) _|_ In2^ 


(23) 


where A 


2 me+2m^ (r—l)(2r—1) 


this theorem is complete. 


and B 


ir^)lZ+2m^) ■ < 2me. The proof of 
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The proof of Theorem [H Clearly, A^g ^ = (1 + rUe + 2m^)N'^ q, and we have known the value 
of A^e.fe of an (r, F)-ugnm N'{t) by Eq.(l5]). For s > 2, the edge number of N'(s) has three 
parts: the number of newly added edges of Y(s), the number A^e,s-i ~Pr^'e,s-i — (1 ~Pr)^e,s-i of 
randomly removed edges of N'{s — 1) with probability Pr and the number pa{l — Pr)-^es-i 
edges randomly added to the remainder obtained from N'(s — 1) by randomly removing edges of 
N'{s — 1) according to probability Pa- Let p = (1 — Pr){l +Pa)- So, we have a recursive formula 

= (24) 

Note that N'(t) and N{t) have the same number of new vertices added and the same number of 
new edges added to — 1) and — respectively. By Eq.© we use repeatedly the recursive 
formula (j24l) . finally, we obtain 
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We compute the edge-cumulative distribution of N'{t) for 0 < r < t as follows: 
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where t = t - h{r)lnk and h{r) = + in ( 2r + m .„) - Clearly, Pe-cum(^) obeys the power law 

distribution and is independent of p = (1 — pr)(l -|-pa)- 
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Notice that Pe-cum{k) is equivalent to Pe-cumik) by Ea. (fT5]) and Eg. ([261) . also -Pe-cum(^) oa 
2rPg_c^m(^)- Because both models N{t) and N'{t) have the same growth operation and the same 
selection of bound-edges which leads to there are similar giant components in them. So, combining 
with (fT6]) gives us 

P(k) 0. 

ZiTTby 

which enables us to conclude that the randomized EBGN-model N'{t) is scale-free. The proof of 
Theorem [1] is complete. □ 




